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We study heavy ion collisions at strong 't Hooft coupling using AdS/CFT correspondence. 
According to the AdS/CFT dictionary heavy ion collisions correspond to gravitational shock 
wave collisions in AdSs. We construct the metric in the forward light cone after the collision 
perturbatively through expansion of Einstein equations in graviton exchanges. We obtain 
an analytic expression for the metric including all-order graviton exchanges with one shock 
wave, while keeping the exchanges with another shock wave at the lowest order. We read 
off the corresponding energy-momentum tensor of the produced medium. Unfortunately 
this energy-momentum tensor does not correspond to ideal hydrodynamics, indicating that 
higher order graviton exchanges are needed to construct the full solution of the problem. 
We also show that shock waves must completely stop almost immediately after the collision 
in AdSs, which, on the field theory side, corresponds to complete nuclear stopping due 
to strong coupling effects, likely leading to Landau hydrodynamics. Finally, we perform 
trapped surface analysis of the shock wave collisions demonstrating that a bulk black hole, 
corresponding to ideal hydrodynamics on the boundary, has to be created in such collisions, 
thus constructing a proof of thermalization in heavy ion collisions at strong coupling. 

Subject Index: 228,121 

§1. General Setup: Expansion in Graviton Exchanges 

AdS/CFT correspondence conjectures that the dynamics of M = 4 SU(N C ) SYM 
theory in four space-time dimensions is dual to the type IIB superstring theory on 
AdSsxSsP In the limit of large number of colors N c and large 't Hooft coupling 
A = g 2 N c (with g the gauge coupling constant) such that N c 3> A S> 1, AdS/CFT 
correspondence reduced to the gauge-gravity duality: J\f = 4 SU(-/V C ) SYM theory at 
N c S> A S> 1 is dual to (weakly coupled) classical supergravity in AdS5. Hence the 
gauge theory dynamics at strong coupling, which includes all-orders quantum effects, 
is equivalent to the classical dynamics of supergravity. Instead of summing infinite 
classes of Feynman diagrams in the gauge theory or using other non-perturbative 
methods, one can simply study classical supergravity in 5 dimensions. For a review 
of AdS/CFT correspondence see.® 

Our goal is to describe the isotropization (and thermalization) of the medium 
created in heavy ion collisions assuming that the medium is strongly coupled and 
using AdS/CFT correspondence to study its dynamics. We want to construct a 
metric in AdSs which is dual to an ultrarelativistic heavy ion collision as pictured 
in Fig. [TJ Throughout the discussion we will use Bjorken approximation of the 
nuclei having an infinite transverse extent and being homogeneous (on the average) 
in the transverse direction, such that nothing in our problem would depend on the 
transverse coordinates x 1 , x 2 P 
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Fig. 1. The space-time picture of the ultrarelativistic heavy ion collision in the center-of-mass frame. 
The collision axis is labeled x 3 , the time is x°. 



We start with a metric for a single shock wave moving along a light cone: 
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ds 2 = ^ 
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Here = g °^? 3 , z is the coordinate describing the 5th dimension such that the 
boundary of the AdS space is at z = 0, and L is the radius of S5. According to 
holographic renormalizationP (T__(af~)) is the expectation value of the energy- 
momentum tensor for a single ultrarelativistic nucleus moving along the light-cone 
in the x + -direction in the gauge theory. We assume that the nucleus is made out 

of nucleons consisting of N 2 "valence gluons" each, such that (T (x~)) oc N 2 , and 

the metric (jl-ip has no A^-suppressed terms in it. 

The metric in Eq. (|l-ip is an exact solution of Einstein equations in AdSs: 
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It can also be represented perturbatively as a single graviton exchange between the 
source nucleus at the AdS boundary and the location in the bulk where we measure 
the metric/graviton field. This is shown in Fig. [21 where the solid line represents the 
nucleus and the wavy line is the graviton propagator. Incidentally a single graviton 
exchange, while being a first-order perturbation of the empty AdS space, is also an 
exact solution of Einstein equations. This means higher order tree-level graviton 
diagrams are zero (cf. classical gluon field of a single nucleus in covariant gauge in 
the Color Glass Condensate (CGC) formalism™ ) . 

Now let us try to find the geometry dual to a collision of two shock waves 
with the metrics like that in Eq. (jl- ip . Defining t\{x~) = 2 Sr (T\ (x~ )) and 

t2( x+ ) = ^jhr {T2++{x + )) for the energy-momentum tensors of the two shock waves 

c 

in the boundary theory, we write the metric resulting from such a collision as 
1? 



ds z = ^ 
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dx 2 +t 2 {x^ 
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+ higher order graviton exchanges > (1-3) 
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Fig. 2. A representation of the metric as a graviton (wavy line) exchange between the nucleus 
at the boundary of AdS space (the solid line) and the point in the bulk where the metric is 
measured (denoted by a cross) . 



The metric of Eq. (|1 -3j) is illustrated in Fig. [3l The first two terms in Fig. [3] (diagrams 
A and B) correspond to one-graviton exchanges which constitute the individual met- 
rics of each of the nuclei, as shown in Eq. (jl-ip . Our goal below is to calculate higher- 
order corrections to these terms, which are illustrated by the diagram C in Fig. [3] and 
by the ellipsis following it. Fig. [3] illustrates that construction of dual geometry to a 
shock wave collision in AdSs consists of summing up all tree-level graviton exchange 
diagrams, similar diagrammatically to the classical gluon field formed by heavy ion 
collisions in CGC.®® While classical gluon fields lead to free-streaming final state,® 
as we will argue below, their AdS graviton "dual" will lead to an ideal hydrodynamic 
final state for the gauge theory similar to the one found in,® though at the same 
time different fromP due to non-trivial rapidity dependence in the case at hand. 




_j_ j_ 

~ N 2 <^ ++ > ~ N 2 <%++> 

ABC 



Fig. 3. Diagrammatic representation of the metric in Eq. (|l-3[) . Wavy lines are graviton propagators 
between the boundary of the AdS space and the bulk. Graphs A and B correspond to the metrics 
of the first and the second nucleus correspondingly. Diagram C is an example of the higher order 
graviton exchange corrections. 



§2. Perturbative Solution of Einstein Equations 

We begin by parametrizing the metric as 
L 2 ( 

ds 2 = — < — [2 + g(x + , x~ , z)] dx + dx~ + [ti(x~ ) z 4 + f (x + , x~ , z)] dx~ 2 
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t 2 (x + ) z A + f(x + ,x~,z) dx +2 + [l + h(x + ,x~,z)} dx\ 



dz' 



(2-1) 



where the functions /, /, g, and h are zero before the collision, i.e., for either x + < 
and/or x~ < 0. The exact Einstein equations (|l-2p for /, /, g, and h are rather 
complicated and are not going to be presented here. Instead we are going to solve 
Einstein equations perturbatively. 

To be more specific let us consider in the boundary theory a collision of two 
ultrarelativistic nuclei with large light-cone momenta per nucleon pf,p 2 , and atomic 
numbers Ai and A 2 . Writing 



ti(x )=mS(x ), t 2 (x + ) = n 2 S(x + ) 



(2-2) 



we want to expand Einstein equations in the powers of /ii and fi 2 ^^^^ The two 
scales /ii and fi 2 can be expressed terms of physical parameters in the problem^ 



fi^p+AjA 1 ^, ^ 2 ^p-A 2 2 Al /3 . (2-3) 

Ai and A 2 are the typical transverse momentum scales describing the two nucleip^ 1 
similar to the saturation scales. Note that and \x 2 are independent of N c . As 
follows from a simple dimensional analysis, combined with Lorentz-transformation 
properties of the relevant quantities, the expansion parameters in four dimensions 
would be 



/ii (x ) 2 x + , and [i 2 (x + ) 2 x , 



(2-4) 



such that the expansion is valid only at early times, when these parameters are small. 

Linearizing Einstein equations (|1-2|) in /, /, g, and h we solve the obtained 
system of differential equation to obtain^ 1 

h(x + ,x',z) = h (x + , X -)z 4 + h 1 (x + ,x-)z 6 (2-5) 

where ho and h\ are determined by the causal solutions of the following equations 

(d + d^) 2 h (x + ,x-) = 8t 1 (x-)t 2 (x + ), (2-6) 



d+ d- hi(x + , x ) 



-t l (x-)t 2 (x+). 



(2-7) 



/, /, and g are easily expressed in terms of h(x + ,x~,z) from Eq. f)2-5[) (see^). 

This lowest-order perturbative solution leads to the energy density of the pro- 
duced medium (in the center-of-mass fv&mej^^} 
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where r 



x+ x . The corresponding center-of-mass energy-momentum tensor is 
/ e(r) \ 
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in terms of the x°, x , x 2 , x 3 components. One can see that the longitudinal pressure 
component in Eq. (|2-9p is large and negative. Under boosts the T 00 and T 33 com- 
ponents of the energy-momentum tensor mix. This implies that there is a frame in 
which the energy density is negative Tq < 0. At this point it is not clear whether 
this result presents a problem, as there may be nothing wrong with energy den- 
sity becoming negative for a short period of timeP^ 1 As we will see below, further 
evolution of the energy density with time leads to disappearance of this negative 
energy-density problem. 

Higher-order perturbative corrections to the energy-momentum tensor (|2-9p in 
powers of \i\ and [12 have been found irPS-HD U p to the fourth order in ^ (i.e., up 
to 0(pf/fc),0(M?/i!),0(/ii^)). 

§3. All-Order Resummation in /X2 

The exact solution of Einstein equations for the collision of two shock waves 
involves resummation of both parameters in Eq. (|2-4p to all orders: such calculation 
appears to be very hard to do. Instead one can resum all orders of 112 (x + ) 2 x~ while 
keeping at the lowest order in fi\ (x~) 2 x + . The corresponding diagram is shown 
in Fig. U} in it one resums multiple graviton exchanges with one shock wave (£2), 
while exchanging only one graviton with the other shock wave (ii). By analogy 
with perturbative QCD calculations we will refer to these class of diagrams as to the 
"proton-nucleus" scattering, with the shock wave t\ being the proton and the shock 
wave ti the nucleus. The applicability region of such approximation is defined by 

Hi {x~) 2 x + <C 1, ^ 2 (x + ) 2 x~~l, (3-1) 

which shows that the resummation is applicable to nucleus-nucleus collisions, only 
in the part of the forward light-cone defined by Eq. (|3-ip . 

ti 




Ax*, Ax + 2 Ax*, 



Fig. 4. A diagram contributing to the metric of an asymmetric collision of two shock waves. 
Such resummation was performed iEpJ using the eikonal approximation. The 
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result for the expectation value of the energy-momentum tensor reads^ 1 

N 2 4 / m / u 2 (x + ) 2 #(x+)#(2r) 



(T- 
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N 2 
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2tt 2 [1 + 8 / u 2 (x+) 2 x-] 3/2 ' 
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(3-2c) 
(3-2d) 



Provided the complexity of the problem at hand, the resulting formulas (|3-2p for the 
energy-momentum tensor are remarkably simple! 

Now we can ask a question: what kind of medium is produced in these strongly 
coupled proton- nucleus collisions? Is it described by ideal hydrodynamics, just like 
Bjorken hydrodynamics was obtained in®? In our case the produced matter distri- 
bution is obviously rapidity-dependent, so it is slightly more tricky to check whether 
Eqs. (|3-2p constitute an ideal hydrodynamics, i.e., whether it can be written as 



= {e + p)u»u v 



prf 



(3-3) 



with the positive energy density e and pressure p. rf" v is the metric of the four- 
dimensional Minkowski space-time and is the fluid 4-velocity. 

For the particular case at hand it is easy to see that the energy-momentum 
tensor in Eq. (|3-2p can not be cast in the ideal hydrodynamics form of ()3-3|) . In the 
case of ideal hydrodynamics one has 



{e + p){u + ) 2 > 0. 



(3-4) 



At the same time {T ++ ) in Eq. ()3-2ap is negative definite. Therefore the ideal 
hydrodynamic description is not achieved in the proton-nucleus collisions. We be- 
lieve this result is due to limitations of this proton-nucleus approximation. As we 
will show below, thermalization (black hole production) is inevitable in the collision 
of the two shock waves considered here. Our conclusion is then that thermaliza- 
tion/isotropization of the medium does not happen in the space-time region defined 
by the bounds in Eq. (|3-1|) . What we found in Eq. (|3-2|) is a medium at some in- 
termediate stage, on its way to thermalization at a later time. It appears that one 
needs to solve the full nucleus-nucleus scattering problem to all orders in both \x\ 
and to obtain a medium described by ideal hydrodynamics. 
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§4. Stopping of Nuclei After Collision 

To better understand dynamics of the shock wave collisions let us follow one of 
the shock waves after the interaction. First we "smear" the delta-function profile of 
that shock wave: 



u x -) = ^-e{x-)e{ ai -x~). 

at 



(4-1) 



A A 1/3 

Here ai oc fii 4 oc -4r-, where the nucleus of radius R\ has A\ nucleons in it. The 

Pi Pi 

"++" component of the energy-momentum tensor of a shock wave after the collision 
at x- = at/2 i^P 



{T + + (x + ,x- = at/2)) 



Nl ixt 



2tt 2 at 



(4-2) 



The first term on the right of Eq. (j4-2[) is due to the original shock wave while the 
second term describes energy loss due to graviton emission. Eq. (|4-2p shows that 
(T ++ ) of a nucleus becomes zero at light-cone times (as ~ P2 in the center-of- 
mass frame) 

1 1 



•^stop 



(4-3) 



VJ^ a 2 a\ /g aI /6 ' 

Zero (T ++ ) would mean stopping of the shock wave and the corresponding nucleus. 
The result can be better understood by doing all-order resummation of graviton 
exchanges with one shock wave performed above for "proton-nucleus collisions" P2J 
The full result for the proton's "++" component of the energy-momentum tensor is 



2vr 2 at + (aF 



for < x < at- 



(44) 



Eq. ()4-4p is illustrated in Fig. [5j in which one can see that the proton loses all of its 
light cone momentum over a rather short time. 




Fig. 5. T ++ component of the proton's energy-momentum tensor after the collision as a function 
of light cone time x + (arbitrary units). 



We thus conclude that the collision of two nuclei at strong coupling leads to 
a necessary stopping of the two nuclei shortly after the collision. If the nuclei 
stop completely in the collision, the strong interactions between them are almost 
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certain to thermalize the system, probably leading to Landau hydrodynamics.^ 
(Rapidity-independent Bjorken hydrodynamics® seems to be unlikely after stop- 
ping. Even before stopping the energy-momentum tensor in Eq. (|3-2p is strongly 
rapidity-dependent.) 

§5. Thermalization: Trapped Surface Analysis 

While the exact solution of Einstein equations for the colliding shock waves 
remains elusive, one can infer whether a black hole will be created in the bulk 
following such collisions by performing a trapped surface analysis.^ A trapped 
surface analysis for shock waves with sources in the bulk has been carried out before 
inP^^SJ However, with the resulting trapped surface being centered around the 
sources, it is not clear to what extend the trapped surface is the property of these 
bulk sources, and what happens to the trapped surface when there are no bulk 
sources, as is the case for our shock wave (|l-ip . 

Performing a trapped surface analysis for a shock wave without sources (jl-lj) does 
not allow one to uniquely fix the position of the trapped surface. 2 -^* We therefore 
start with a shock wave having an extended source in the bulk with the only non-zero 
component of the bulk energy-momentum tensor being J = S(x~) 5{z — zq). 

The corresponding metric is^ 1 



L 2 

ds 2 = \-2dx + dx~ +cj)(z)5(x-)dx~ 2 + dx\ + dz 2 } (5-1) 

z ' 



with 



2n 2 EL 2 
N 2 



zz, z <z 



1, Z > Z . 



(5-2) 



As one can readily see, the metric (|5-ip reduces to that in Eq. (|l-ip if we send the 
source to the IR infinity in the bulk by taking zq — > oo limit of the metric in Eqs. 
(|5-ip and (|5-2p keeping fi defined by 

2tt 2 EL 2 , N 

"=ivf-^ (5 ' 3) 

fixedP 

Marginal trapped surface for a collision of two shock waves given by Eq. (|5-ip 
was found inP^ In the zq — > oo, [i fixed limit that trapped surface reduces to 



+ n - V 
x = U, x = — — 



z 4 



M ~ 4 / 3 2- 2 / 3 l (5-4) 



with an analogous expression for the other shock wave obtained by interchanging 
x + -f-7- x~ in Eq. (|5-4p . (For simplicity we work in the center-of-mass frame where 
/^l = = M-) The trapped surface in Eq. (|5-4|) is independent of the shape of the 
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Fig. 6. An illustration of the trapped surface in the collision of two sourceless shock waves. Vertical 
axis is the bulk ^-direction, the horizontal left-right axis can be thought of either as the collision 
axis or as the time direction. The trapped surface is shaded. 



source being sent to deep IR, as shown inP^ The trapped surface from Eq. (|5-4p is 
depicted in Fig. [6j 

The existence of trapped surface proves that gravitational collapse is inevitable. 
That is a black hole will be created in a bulk following a collision of two sourceless 
shock waves. In the boundary theory this means that a thermal medium is created, 
which is described by ideal hydrodynamics. 

The (lower bound for the) produced entropy per unit transverse area A± can be 
found by calculating the area of the trapped surface, which yields^ 

■f - = -\ (2 W ^ 2 ) 1/3 . (5-5) 
A± 2 ir z 

Since the trapped surface analysis does not "know" anything about shock wave 
thickness (e.g. a\), we conclude that the thermalization time is only a function of 
fj,\ and fi2, which gives 

in agreement with the thermalization time suggested inP^ While numerically this 
thermalization time is too short to be relevant for RHIC data, it is parametrically 
shorter than the stopping time (|4-3p . making our model somewhat more relevant for 
description of real-life collisions. As \x\ fi2 ~ s with s the center of mass energy of 
the collision, we get 

4- « * 1/3 (5-7) 
A± 

in agreement with the result obtained inP^ Finally note that, since, as follows 
from Eq. ()2-3p . /zi and \ii are iV c -independent, the produced entropy scales oc N^, in 
agreement with iV c -counting in a perturbative QCD calculation of particle production 
for a collision of two nuclei with "valence gluons" in their nucleons. 
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